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Using recent measurements of the spectrum and chemical composition of the highest energy 
cosmic rays, we consider the sources of these particles. We find that the data strongly prefers 
models in which the sources of the ultra-high energy cosmic rays inject predominantly intermediate 
mass nuclei, with comparatively few protons or heavy nuclei, such as iron or silicon. If the number 
density of sources per comoving volume does not evolve with redshift, the injected spectrum must be 
very hard (a ~ 1) in order to fit the spectrum observed at Earth. Such a hard spectral index would 
be surprising and difficult to accommodate theoretically. In contrast, much softer spectral indices, 
consistent with the predictions of Fermi acceleration (q ~ 2), are favored in models with negative 
source evolution. With this theoretical bias, these observations thus favor models in which the 
sources of the highest energy cosmic rays are preferentially located within the low-redshift universe. 


I. INTRODUCTION 

Despite considerable advances in the measurement of the spectra, arrival directions, and chemical composition of the 
ultra-high energy cosmic rays (UHECRs), we remain ignorant of the sources of these particles. In particular, although 
it had been hoped that anisotropies in the arrival directions of the UHECRs would be detectable by the Pierre Auger 
Observatory (PAO), such studies have thus far not conclusively identified a correlation between these particles and 
any known class of astrophysical objects [T]. Efforts to identify UHECR sources using secondary signatures have also 
been inconclusive [a El- No signal capable of revealing the origin of the UHECRs has yet appeared. 

Despite the lack of such a definitive result, the available information bearing on this puzzle has significantly increased 
over the past few years, allowing for considerable progress to be made. Of particular importance are the PAO’s 
measurements of the chemical composition of the UHECR spectrum HE], which reveal the highest energy cosmic 
rays to be dominated by intermediate mass nuclei, as opposed to either protons or heavy nuclei, such as iron or silicon. 
While earlier data also supported similar conclusions [5], the PAO’s most recent template-based composition study [S] 
has made a rather detailed determination possible. 

Considering the most recent spectrum and composition measurements from the PAO H above an energy threshold 
of 10^® ® eV, we investigate models for the injected spectrum, composition, and redshift distribution of the UHECR 
sources. Our main finding is that in models with no significant evolution with redshift (a constant number of sources 
per comoving volume), the injected spectrum from the sources of the UHECRs must exhibit a very hard spectral index, 
a ~ 1.1. Such a hard spectrum would be surprising and difficult to accommodate theoretically (although scenarios 
have been proposed mu)- Models with positive redshift evolution (a greater number of sources per comoving volume 
at high redshift) only exacerbate this problem. In models with significant negative evolution, however, the measured 
spectrum can be accommodated for more well-motivated values of the injected spectral index, near that generically 
predicted from Fermi acceleration, a ~ 2. This theoretical bias for soft cosmic ray injection spectra supports the 
conclusion that the sources of the UHECRs are predominantly located in the low-redshift universe. One notable 
implication of this result is that a low cosmogenic neutrino flux is predicted, well below the expected sensitivity of 
existing or planned experiments. 

In this paper, we investigate the propagation of ultra-high energy cosmic ray nuclei and compare the predictions of 
various models to the current data from the PAO. In Sec.|^ we describe our treatment of UHECR nuclei propagation 
and the range of models considered. In Sec. m we describe our results, which favor models with negative redshift 
evolution and which inject significant quantities of intermediate mass nuclei, such as helium, carbon, nitrogen, and 
oxygen. In Sec. |IV[ we discuss the contribution from UHECR propagation to the diffuse gamma-ray background and 
the flux of cosmogenic neutrinos. Lastly, in Sec. |V] we summarize our results and conclude. 
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II. THE PROPAGATION OF ULTRA-HIGH ENERGY COSMIC RAY NUCLEI 

To test various models for the sources of the highest energy cosmic rays, we employ a Monte Carlo description 
of UHECR propagation, as originally described in Ref. [12] ■ In this calculation, UHECR protons and nuclei are 
propagated through the cosmic microwave background (CMB) and the cosmic infrared background (CIB), undergoing 
photo-disintegration, photo-pion interactions, pair production, and redshift energy losses. In order to allow for a large 
number of models to be studied quickly, we inject particles with energies and distances over discretized logarithmic 
ranges, and then sum these results with appropriate weighting factors, so as to generate results for more general 
functions of the source energy spectra and spatial distribution. 

Eor each model, we calculate the average depth of shower maximum, (Amax), and its RMS variation, RMS(Ai„ax), 
as a function of energy, and compare this to the values measured by the PAO |4]. In order to encapsulate the 
uncertainty associated with the modeling of the shower development, we plot a band which encompasses the values 
as determined using the different hadronic models, EPOS-LHC [T3] and QGSJET-II-04 [HI, which have each been 
updated to incorporate the first results from the Large Hadron Collider, as well as SIBYLL m- The cross sections 
and target photon spectral distributions relevant for proton energy losses are sufficiently well understood so as not 
to introduce large uncertainties in our results. Eurthermore, although the uncertainties regarding our knowledge 
of the photo-disintegration cross sections and the CIB spectral distribution relevant for nuclei propagation are not 
necessarily negligible, such uncertainties are unlikely to qualitatively impact our conclusions. In our calculations, we 
adopt the cross sections as given in Ref. [T6| and the CIB spectral distribution as described in Ref. m- 

In this study, we consider models with a wide range of redshift distributions, energy spectra, and chemical compo¬ 
sitions. Civen our present ignorance of the nature of these sources, we consider a simple source evolution model with 
the aim to encapsulate a broad set of possible histories, which we parameterize by: 

OC (1 -I- z)”, z < Zmax, (1) 

aVc 

where dN/dVc is the number of sources per comoving volume. We will present results for a wide range of evolution 
models, from n = 3 to n = —6. We adopt a maximum redshift value of Zmax = 3, which is well beyond the distance 
that UHECRs are able to propagate. 

Eor the spectra of UHECRs injected from sources, we adopt the following form: 

diV 

— ocU “ exp[-E;/E;z.niax], (2) 

where a describes the spectral index and i?z,max = iZj26) x Upe.max is a rigidity dependent maximum energy. 

Eor the chemical composition of the UHECR spectrum at injection, we consider an arbitrary mixture of 5 nuclear 
species, namely: protons, helium, nitrogen, silicon, and iron. Physically, we take these categories to represent groups 
of nuclei. For example, the nitrogen fraction can be thought of as a proxy for the fraction of injected cosmic rays that 
are carbon, nitrogen, or oxygen. 


HI. FITS TO THE MEASURED SPECTRUM AND CHEMICAL COMPOSITION 

In this section, we describe the statistical methods and main results of our analysis, including the range of UHECR 
models that provide a good fit to the combined spectrum and composition measurements from the PAO. Given 
the large number of free parameters being varied in our analysis, we adopt a method that makes use of a Markov 
Chain Monte Carlo (MCMC) routine, employing the Metropolis-Hastings algorithm to scan through the spectral 
and composition likelihood landscape and determine the regions that provide a good fit to the data. This method 
allows local minima regions within the parameter space to be explored and “escaped from” such that the true global 
minimum can be effectively located. In these scans, the Likelihood function (for given values of n and Zmax) is given 
by: 

j /hcj /n ; /sij L^Fe,max; Q:) OC exp( y (/p; /He; /n ; /si; A’Pe.max; Q;)/2) , 

where the quantities, fi, denote the fraction of the injected cosmic ray spectrum that consists of the nuclear species, 
i. Note that the parameter /pe is not free, but is instead fixed by the constraint, /p -I- /ue + /n + /si + fpe = 1- We 
only include data at energies above 10^® ® eV in our fit, as it is not clear that the cosmic ray spectrum is dominated 
by extragalactic particles at lower energies [HI [H] . Systematic errors in the analysis have been added in quadrature 
to statistical errors. Throughout our analysis, we adopt the following systematic errors: 14% on energy, 10 g cm“^ 
on (Ainax), and 2 g cm”^ on RMS(Ainax) [Ij- 
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FIG. 1: The best-fit models for a source evolution model dN/dVc oc (1 -I- a)", up to 2max = 3 with different indices: n = 3, 
n = 0, n = —3 and n = —6 from top to bottom. In the left frame, we compare the total predicted UHECR spectrum to 
that measured by the Pierre Auger Observatory. The dashed and dotted curves in this frame denote the contribution from 
individual nuclear mass groups. In the middle and right frames, we compare the prediction of this model to the depth of shower 
maximum (Xmax) and its RMS variation, again as measured by Auger [3]. 


Figurej^shows the best-fit spectra, (Xmax) and RMS(Xmax) distributions for four different evolution indicies: n = 3, 
0, —3, and —6. We begin by considering the case with no source evolution (n = 0). In this case, the best-fit is found 
for sources which inject a large fraction of helium (/hb = 0.53) and particles in the nitrogen group (/n = 0.29), along 
with smaller but non-negligible quantities of protons {fp = 0.17) and very few heavier nuclei (/si = 0.0, fpe = 0.01). 
As mentioned in the introduction, the injected spectrum is very hard in this model, with an index of a = 1.1 and 
^'Fe.max = 10^°'^ eV. This is in considerable contrast to the softer spectra generally predicted by Fermi acceleration, 
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FIG. 2: Frequency plots from the MCMC scan of the spectral parameters a and i?Fe,max for four values of the source evolution 
parameter, n. Models with no evolution (n = 0) or positive evolution (n = 3) favor very hard spectral indices for the injected 
spectrum, whereas models with negative evolution (n = —3 & n = —6) prefer significantly softer injected spectra. 


a ~ 2. The results for the best-fit model with no source evolution are shown in Fig. 

We also find good fits for positive (n = 3) and negative (n = —3 & —6) source evolution as shown in the other 
panels of Fig. in each case favoring models with large fractions of helium and nitrogen at injection and ifFe.max — 
(1.5 — 3) X 10^ eV (see Table]^. The main difference among the best-fit models is the value of the injected spectral 
index, which varies from a =0.6 for u = 3 to a =1.8 for n = —6. In Fig. we compare the best-fit models for 
the four different evolution models to the mass composition inferred by PAO via the template-based fit of the Xmax 
distribution The vertical dashed line indicates the low energy threshold of our fit. For all evolution models 
the observed best-fit mass spectrum after propagation agrees with the PAO observation within uncertainties. In 
particular, heavy mass groups beyond nitrogen have only a small contribution in the observed spectrum. 

Moving beyond best-fit models, we show in Fig. [^the weighted distributions of the spectral parameters a and ifpe.max 
that are found by our MCMC scan, for each evolutionary model. Although a modest shift in Ape.max is observed with 
n, the greatest variation is observed among the values of a that are favored. Models with no evolution or positive 
evolution favor very hard spectral indices for the injected spectrum, whereas models with negative evolution prefer 
softer injected spectra. In Table we also provide the posterior mean and standard deviation for these parameters, 
as well as for the injected composition fractions. 

In addition to the spectral parameterization described in Eq. we have also considered spectra with sharper 
than exponential cutoffs. Such spectra are motivated from particle acceleration scenarios in which electrons are both 
accelerated and radiatively cooled within the source environment Eniin]- No qualitative changes to our conclusions 
were found in this case, however. 

It should be noted that the distribution of the spectral parameters presented in Fig. are correlated, and are 
thus not independently constrained. As pointed out previously in Ref. [22] (which considered the specific case of 
n = 3), the models with the largest values of a also exhibit large values of Epe^max- Since the time of that study, the 
systematic uncertainties associated with the PAO’s measurements have been reduced considerably, and updates from 
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n = —6 


n = —3 


n = 0 


n = 3 

Parameter 

Best-fit 

Value 

Posterior Mean & 
Standard Deviation 

Best-fit 

Value 

Posterior Mean & 
Standard Deviation 

Best-fit 

Value 

Posterior Mean 
Standard Deviation 

Best-fit 

Value 

Posterior Mean Sz 
Standard Deviation 

fp 

0.03 

0.14 ±0.12 

0.08 

0.15 ±0.13 

0.17 

0.17±0.16 

0.19 

0.20 ±0.16 

/Ue 

0.50 

0.21 ±0.17 

0.42 

0.17±0.16 

0.53 

0.20 ±0.17 

0.32 

0.23 ±0.20 

/n 

0.40 

0.50 ±0.18 

0.42 

0.51 ±0.19 

0.29 

0.47 ±0.19 

0.43 

0.45 ±0.21 

/si 

0.06 

0.11 ±0.12 

0.08 

0.12 ±0.13 

0.0 

0.11 ±0.12 

0.06 

0.078 ± 0.086 

fPe 

0.01 

0.052 ± 0.039 

0.0 

0.053 ±0.042 

0.01 

0.050 ± 0.038 

0.0 

0.044 ±0.034 

Q 

1.8 

1.83 ±0.31 

1.6 

1.67 ±0.36 

1.1 

1.33 ±0.41 

0.6 

0.64 ±0.44 

logio(%^) 

20.5 

20.55 ±0.26 

20.5 

20.52 ±0.27 

20.2 

20.38 ±0.25 

20.2 

20.16 ±0.18 


TABLE I: The parameters which provide the best-fit to the spectrum, (Xmax), and RMS(Xmax), as measured by the Pierre 
Auger Observatory, for four choices of the source evolution model. Models with no evolution (n = 0) or positive evolution 
(n = 3) favor very hard spectral indices for the injected spectrum, whereas models with negative evolution (n = —3, -6) prefer 
significantly softer injected spectra. We also show the posterior mean value and standard variation of the injected spectrum 
and composition parameters resulting from our MCMC scan. 


recent LHC runs have reduced in the variation between the predictions of various hadronic models. After taking these 
improvements into account, it is apparent that the scenarios considered in Ref. [22) . which were previously found to 
be in reasonable agreement with the data, are now disfavored. 

We also note that the proton to helium ratios in Table within their 1 sigma standard deviation regions, appear 
unphysically small. This feature appears to be an artifact of the energy threshold of 10^® ® eV adopted. Indeed, 
a reduction of this energy threshold in the analysis allows much larger ratios. Such a decrease of the threshold 
energy, however, potentially also demands a more complicated setup. With the composition of the flux in the energy 
range below 10^® ® eV and above 10^® eV being indicated to be light, we highlight these complications through the 
consideration of an additional proton flux component, whose spectral index/source evolution is not tied to that of 
the heavier nuclei, to account for this flux. We find that for reasonable fits to the spectrum to be obtained, either 
softer spectral indices than for the > 10^® ® eV component (as found previously by others m), or a stronger source 
evolution are required for this additional light component. 


IV. CONSTRAINTS FROM COSMOGENIC GAMMA-RAYS AND NEUTRINOS 

In the previous section, we compared the predictions of a wide range of UHECR models (scanned by MCMC) to the 
spectrum and composition measurements reported by the Pierre Auger Collaboration. In some respects, these results 
were found to be largely insensitive to variations in the source evolution parameter, n. In particular, we find good fits 
to the combined spectrum and composition data for each evolution scenario, and the injected chemical composition 
preferred by the fit is only mildly impacted by our choice of the evolution model. In other ways, however, our results 
depend very strongly on the choice of source evolution. Most notably, the injected spectral index required to fit the 
spectrum measured by the PAO varies considerably with evolution, favoring a < 1 for models with significant positive 
evolution and a ~ 2 for models with significant negative evolution. 

In addition to the impact on the spectral index the flux of cosmogenic gamma-rays and neutrinos from the propaga¬ 
tion of UHECR varies considerably with the evolutionary model being considered. Interactions of the UHECRs with 
the cosmic microwave and infrared backgrounds invariably generate a population of high-energy neutrinos, known 
as the cosmogenic neutrino flux. In particular, for proton-dominated models, photo-pion production with the CMB 
leads to a cutoff in the spectrum at Eqzk — 50 EeV [23l - f25) . The resulting neutrino flux becomes maximal at a few 
EeV [? ]. 

If the UHECR spectrum is dominated by intermediate mass nuclei, as indicated by cosmic ray observations by the 
Auger Collaboration, however, the resonant interaction of cosmic ray nucleons with the CMB is shifted to higher 
energies, (A/56) x 3 ZeV, where A is the atomic mass. For heavy or intermediate mass nuclei, the energy at which 
neutrinos are most efficiently generated is shifted upward, to where significantly fewer cosmic rays exist. As a result, 
nuclei-dominated cosmic ray models predict significantly lower fluxes of cosmogenic neutrinos than proton-dominated 
scenarios. 

Due to this increased threshold of cosmogenic neutrino production for cosmic ray nuclei, cosmic radiation at infrared. 
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FIG. 3: The predicted chemical composition of the best-fit UHECR spectrum at Earth compared to that reported by the Auger 
Collaboration [^. Results are shown for a source evolution model dN/dVc oc (1 -I- z)'^, up to Zmax = 3 with different indices: 
n = 3, n — 0, n = —3 and n = —6 from top to bottom. For these results the EPOS-LHC [13] hadronic model is adopted and 
the data and model line colors indicate the nuclear species range applicable. 


optical, and ultra-violet wavelengths can be important targets, and are included in most modern cosmogenic neutrino 
calculations [I1IM33. In general, the peak of the cosmogenic neutrino spectrum predicted for nuclei-dominated 
models is shifted downward, to the 1-100 PeV range, and the overall normalization is reduced relative to proton- 
dominated models, generally below present experimental sensitivities. The overall cosmogenic neutrino flux also 
depends on the maximal energies injected and the redshift evolution of UHECR sources. An estimate of the lower 
limit of these pessimistic models was given in Ref. [51] . 

In the right panel of Fig. we show the cosmogenic neutrino flux (summed over all flavors) for the best-fit 
models for four choices of redshift evolution (n = 3, 0, —3, —6). The higher energy peak at EeV energies is due 
to neutrinos produced on the CMB photons, whereas the lower peak at PeV energies is largely due to interactions 
with the infrared background m- These predictions are well below the present limits placed by the IceCube [ 55 ] . 
Anita [ 55 ] and PAO [ 57 ] Collaborations, and are even too low to be reached by future radio Cherenkov observatories, 
such as ARA [ 55 ], or ARIANNA [ 55 ]. For orientation we show the projected ARA-37 3yr sensitivity in the plot. 
These predicted cosmogenic neutrino fluxes are also much lower than the extraterrestrial neutrino flux reported by 
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the IceCube Collaboration, which is near the level of E^Jl^ ~ (9.5 ± 3.0) eV cm^ sr“^ over the energy range of 
~10 TeV-PeV [40]. It has been speculated that this emission could be indirectly related to the flux of UHECRs via 
calorimetric processes in cosmic ray reservoirs affecting the low energy tail of the emission spectrum mi 132] . Such 
environments have been discussed within the context of starburst galaxies US] |33], Galactic outflows |3SJ US], or 
clusters of galaxies [37ti3^ . 

The photo-hadronic interactions with the cosmic radiation backgrounds responsible for the flux of cosmogenic 
neutrinos also generate a flux of high-energy gamma-rays. In addition, non-resonant Bethe-Heitler pair production 
via cosmic ray scattering off CMB photons creates a population of high energy electrons and positrons. However, 
these electro-magnetic contributions undergo rapid cascades via pair production and inverse-Compton scattering in 
the CMB, resulting in a cascaded sub-TeV gamma-ray spectrum. These gamma-ray contributions are shown for the 
best-fit models in the left plot of Fig.j^in comparison to the isotropic gamma-ray background inferred by Fermi [SO]. 

For the best-fit model with no source evolution (n = 0, Zmax = 3), the electromagnetic cascades resulting from 
UHFCRs produce approximately 10% of the measured isotropic gamma-ray background above 10^^ eV [SO]. In models 
with positive evolution (n = 3, Zmax = 3), the contribution increases to 20%, while models with negative evolution 
contribute significantly less to the diffuse gamma-ray flux. This result is not surprising, as a reduction in the evolution 
parameter amounts to moving the sources of the UHFCRs to more local distances, reducing their mean propagation 
time through the extragalactic radiation fields that give rise to these losses (predominantly the CMB). The entirety of 
the isotropic gamma-ray background observed by Fermi can be easily accommodated by the sum of known gamma-ray 
source classes (blazars, starforming galaxies, and radio galaxies) [5TH55] . and Fermts most recent measurement leaves 
relatively little room for an additional component originating from UHECR propagation (see, for example. Fig. 3 of 
Ref. [51]). However, a contribution at the level of ~10-20% at energies above 10^^ eV is likely compatible with the 
current Fermi data. 

We also briefly consider the cascade contribution from a separate dominant extragalactic proton component in 
the range ® eV for different source evolution models. We find that such a setup gives rise to a factor of 

~2-4 increase in the cascade flux levels, over our previous results, for the different evolution scenarios (n = 3, —6) 
considered. The increase by a factor of ~ 2 was found for the n = —6 evolution and the increase by a factor of 
^ 4 for the n = 3 evolution, with intermediate evolution scenarios receiving intermediate flux increases. Though 
such increases still sit below the isotropic background level, this result highlights the importance in determining the 
transition energy at which the extragalactic sources start to dominate the cosmic ray spectrum, as already highlighted 
by others [Ml liS] . 

We also note that the model-to-model variation in the contribution to the isotropic gamma-ray background predicted 
by our MCMC scan is quite small (when the source evolution model is fixed). For the four source evolution models 
considered, we find that the predicted contribution to the gamma-ray background is always less than 20%. 


V. SUMMARY AND DISCUSSION 

In this study, we have used a Markov Chain Monte Carlo approach to explore a wide range of models for the 
sources of the highest energy cosmic rays, varying the injected spectra, chemical composition, and redshift distribution. 
After comparing the predictions of these models to the most recent data from the Pierre Auger Observatory [3], we 
find that models dominated by intermediate mass nuclei are strongly favored. This is consistent with a template- 
based composition analysis of the ATmax distribution [^. Among these models, those without redshift evolution 
(corresponding to a constant number of sources per comoving volume, n = 0) are capable of fitting the observed 
spectral shape only if the injected spectrum is very hard, with a spectral index of a ~ 1.3. If the number of sources 
increases at low redshifts (n < 0), however, softer spectra consistent with the expectations of Fermi acceleration 
(a ~ 2) are favored. 

We argue that these considerations significantly favor UHECR models with negative source evolution, and strongly 
disfavor positive evolution scenarios. These findings both strengthen previous results indicating a local proximity of 
UHECR sources [56], and further suggests a reduction in the source comoving density at higher redshifts. 

We also show that the contribution of electro-magnetic cascades from the production of high-energy gamma-rays, 
electrons and positrons during the propagation of UHE CRs are consistent with the isotropic gamma-ray background 
inferred by the Fermi Gamma-Ray Space Telescope |50j . The flux of cosmogenic neutrinos is also consistent with 
upper limits. In general, a strong negative evolution, i.e. a stronger contributions of local sources decreases the flux 
of these contribution and makes a direct observation more challenging. 

In light of these results, it is interesting to contemplate which classes of UHECR sources could plausibly be dis¬ 
tributed preferentially within the low-redshift universe. One interesting source class are low-luminosity gamma-ray 
BL Lacertae (BL Lac) objects. While Fermi measurements have revealed that the number density of bright BL Lacs 
peaks at a fairly high redshift of z ~ 1.2, the more numerous low-luminosity erg s“^) and high-synchrotron 




FIG. 4: The contribution of cosmogenic gamma-rays (left plot) and cosmogenic neutrinos (all flavor; right plot) from UHECR 
propagation for the best-fit models in four different source evolution scenarios (n = 3, 0, —3, —6). For comparison, we also show 
the contribution resulting from a purely-proton model with no source evolution (a = 2.05, i?max,p = 10^®'® eV, n = 0). The left 
panel shows the isotropic gamma-ray background inferred by Fermi [50] . All best-fit models considered here are consistent with 
the gamma-ray background. In the models favored by the data (nuclei-dominated, n ^ 0), also very low cosmogenic neutrino 
fluxes are predicted. The thick dotted line in the right panel shows the sensitivity of the proposed ARA-37 after three years of 
observation |38| . 


peaked members of this population exhibit negative source evolution, and thus are overwhelmingly distributed at low 
redshifts m- And while many other possibilities remain, the results of this paper provide some support for this class 
of objects as the sources of the UHECRs. Given the modest number of low-redshift BL Lacs observed by Fermi, 
the lack of anisotropy in the UHECR spectrum would likely require significant deflection by magnetic fields, whether 
throughout the intergalactic medium, or within the volume of local group or Milky Way. 

Lastly, the contribution of an additional extragalactic proton component, to explain the flux below 10^® ® eV 
and above 10^® eV, was also considered. The composition of this component, indicated to be light, motivating the 
consideration for it to be extragalactic in origin. Curiously, however, this component was found to require either a 
softer source spectral index, as also noted by others previously [19], or a stronger source evolution parameter to the 
higher energy component we have focused on. The additional gamma-ray cascade contributions from this component 
was found to be a factor of ~2-4 larger than that from the > 10^® ® eV contributions considered. This result highlights 
the importance in the determination of the transition energy at which extragalactic sources start to dominate the 
cosmic ray flux. 
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